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We study two-dimensional interacting electrons in a weak perpendicular magnetic field with the 
filling factor S> 1 and in the presence of a quenched disorder. In the framework of the Hartree-Fock 
approximation, we obtain the mean-field phase diagram for the partially filled highest Landau level. 
We find that the CDW state can exist if the Landau level broadening 1/2t does not exceed the 
critical value l/2rc = 0.038ujh- Our analysis of weak crystallization corrections to the mean-field 
results shows that these corrections are of the order of <^ 1 and therefore can be neglected. 
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I. INTRODUCTION 

A two-dimensional electron gas (2DEG) in a perpen- 
dicular magnetic field was a subject of intensive studies, 
both theoretical and experimental, for several decades. 
The behaviour of the system in a strong magnetic field 
where only the lowest Landau level is occupied, has been 
investigated in great details^. Several attempts2*^ were 
made in order to incorporate the case with larger fill- 
ing factors V > 1 into the theory. Usually in these ap- 
proaches, the ratio of the characteristic Coulomb energy 
(at distances of the order of the magnetic length) to the 
cyclotron energy, has been assumed to be small. How- 
ever, in a weak magnetic field this is not the case, and the 
characteristic Coulomb energy exceeds the cyclotron en- 
ergy. An attempt to investigate the situation with large 
Coulomb energy was made in Refi4. 

The progress in understanding of the clean 2DEG in a 
weak magnetic field was achieved by Aleiner and Glaz- 
mar;4. They have derived the low-energy effective the- 
ory on the partially filled highest Landau level by using 
the small parameter 1/v ^ 1. By treating the effec- 
tive interaction within the Hartree-Fock approximation, 
Koulakov, Fogler, and Shklovskii^ predicted a unidirec- 
tional charge-density-wave (CDW) state (stripe phase) 
for the half-filled highest Landau level at zero temper- 
ature and in the absence of disorder. Moessner and 
Chalker4 extended the ideas of Fukuyama, Platzmann 
and Anderson'' to the case of a partially filled highest 
Landau level and showed the existence of the mean-field 
CDW state on the half-filled Landau level below some 
temperature Tq- 

Recently, the existence of compressible states near half- 
filling with anisotropic transport properties was demon- 
strated experimentally for high Landau levels^' This 
stimulates an extensive study of the clean 2DEG in a 
weak magnetic field and pinning of stripes by disorder—. 

In the clean case, the properties of the CDW states 
can be described on the basis of the low energy effective 
theory for smooth "elastic" deformationaii. Recently, at- 



tempts were made to derive such a theory microscopically 
starting from the mean-field solutioniSii^. The effects of 
a quenched disorder on the unidirectional CDW state 
(stripe phase) were investigated in the framework of the 
phcnomenological elasticity theoryi^, and a rich variety 
of different regimes, which depend on the strength of dis- 
order, were found. However, to identify the phcnomeno- 
logical parameters of the theory, a successive microscopic 
theory should be developed. 

At present, a thorough microscopic analysis of the ef- 
fects of disorder on the mean-field transition from the 
uniform state to the CDW one, as well as on the phase 
diagram of the mean-field CDW states, is absent. The 
main objective of the present paper is to investigate these 
effects on the existence of the mean-field CDW states in 
2DEG in a weak perpendicular magnetic field H (filling 
factor V 1). For the considered case of a large number 
of the occupied Landau levels, the mean- field analysis is 
legitimate because the fluctuations of the order param- 
eter are strongly suppressed^^. On the other hand, the 
mean-field approach cannot be applied to the critical re- 
gion in the direct vicinity of the phase transition. This 
region is however small and does not lead to any signif- 
icant uncertainty in our results for critical temperatures 
of the transitions. 

We assume the presence in the system of a weak 
quenched disorder, i.e. the elastic collisions time satis- 
fies the condition tq 3> ^ where ujh = eH/m is the 
cyclotron frequency, e the electron charge, and m the ef- 
fective electron mass. Under this condition, the Landau 
level broadening 1/2t, which is of the order of ^JuJhTq/tq, 
is much less than the spacing loh between them. This 
case can be realized in high mobility samples which 
were used for experimental studies of the anisotropic 
magnetotransportSiiflii^. Keeping in mind that the re- 
lation 1/t is usually hold, one expect a much pro- 
nounced influence of the quenched disorder on the prop- 
erties of electrons on the partially filled highest Landau 
level even for a small level broadening 1/2t <C ujh- 

One of the main results of our paper is that at zero 
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temperature the mean-field CDW state is destroyed wlien 
the Landau level broadening exceeds the critical value 
1/2tc — ATq/ti. At nonzero temperatures the quenched 
disorder leads to the decrease of the temperature of the 
CDW instability as compared to the clean case. The 
physical reason is that the scattering on impurities breaks 
the CDW correlations, and therefore results in the de- 
struction of the coherent CDW state. This is somewhat 
similar to the suppression of the critical temperature in 
conventional superconductors by magnetic impuritiesA^ 
or in anisotropic superconductors by nonmagnetic impu- 
rities^. 

The paper is organized as follows. In Sec. we in- 
troduce the formalism that allows us to evaluate the free 
energy of the CDW state in the presence of disorder. In 
Sec, mil we investigate the instability of the uniform state 
towards the formation of the CDW state, and present 
the mean-field phase diagrams at the half-filling and ar- 
bitrary temperature, and at zero temperature and arbi- 
trary filling. The weak crystallization corrections to the 
mean-field solution are presented in Sec. II VI Sec.Mcon- 
tains the comparison of the theory with the recent exper- 
imental and numerical results. We end with conclusions 
in Sec. |Vl| 

II. FREE ENERGY OF THE CDW STATES 

We consider two-dimensional interacting electrons in 
the presence of a weak quenched disorder and a weak 
perpendicular magnetic field. The parameter that char- 
acterizes the strength of the Coulomb interaction is 
Ts = ^/2e'^ /evp with vp being the Fermi velocity and 
e the dielectric constant of a media. We assume that 
the Coulomb interaction between the electrons is weak. 



Ts ^ 1, and the magnetic field obeys the condition 
Nvs ^ 1, where N = [v/2\ is the integer part oi v/2. In 
this case it is possible to construct an effective field the- 
ory for the electrons on the highest partially filled Landau 
level by integrating out all other degrees of freedom^ti^. 
We also assume that the electrons on the partially filled 
highest Landau level are spin-polarized. This assumption 
is based on the calculations'^' that shows the existence 
of fractional states, composite fermions, and skyrmions 
only on the lowest and the first excited Landau levels, as 
well as on the experimental observations. 

In order to study the transition from the uniform state 
to the CDW one we employ the Landau expansion of the 
free energy in powers of the CDW order parameter A(qj), 
where the vectors ^.j that characterize the CDW state, 
have the same length^ qj = Q. We perform the expansion 
up to the forth order in the CDW order parameter under 
the assumption TVr^ 3> 1. In this case the Hartree-Fock 
approximation is well justified® because the corrections 
are small in the parameter ub/Ih — V-^^s ^ 1' where 
qb = e/me^ is the Bohr radius and Ih = \/ y/mujH the 
magnetic length. 

A. Formalism 

The thermodynamical potential of the spin-polarized 
2DEG projected on the iVth Landau level in the presence 
of the random potential Vrfii (r) and the magnetic field is 
given by 

^ = --^ J '^bP^^] J 'D[Vd^s]V[Vd^s] exp(5[V^,^,Fd,,]). 

(1) 

where the action S[ip, ipj Vdis] in the Matsubara represen- 
tation has the form 



5 = 



I UJ„ 



-fj,~Ho~ Vdts (r) 



T 
2" 



E 



(r)C„ - (r)t/o(r, r (r)C„ 



+ ^1 



(r 



(2) 



Here (r) and "0"^ (r) are the annihilation and creation 
operators of an electron on the iVth Landau level, T the 
temperature, /i the chemical potential, w„ = TTT{2n+ 1) 
the Matsubara fermionic frequency, and Vn = 'i.irTn the 
bosonic one. The free Hamiltonian TLq for 2D elec- 
trons with mass m in the perpendicular magnetic field 
H = €abdaAb is Ho = (-iV - eA)'^/{2m). The screened 
electron-electron interaction [/o(r) on the A^th Landau 
level takes into account the effects of interactions with 
electrons on the other levels, and has the form (see Refi^ 
for the clean case and Ref;-- for the weakly disordered 



case) 



Uo{q) = 7^ — 7 K , 

qaB V GujhtJ 

(3) 

where Rc ~ Ih\/v is the cyclotron radius on the A^th 
Landau level and Jo{x) the Bessel function of the first 
kind. The range of the screened electron-electron inter- 
action (j2Jl is determined by the Bohr radius ub ■ We also 
assume the Gaussian distribution for the random poten- 
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tial Vdis{r) 

V[Vd^s{r)] = 



/7rg 



exp 



(4) 



where g = l/npTQ, p is the thermodynamical density of 
states, and introduce Nr replicated copies of the system 
labeled by the replica indices a = l,...,Nr in order to 
average over the disorder. 



B. The Hartree-Fock decoupling and the average 
over disorder 



The CDW ground state is characterized by the order 
parameter A(q) that is related to the electron density 



(p(q)) = L^LynLFNiq)A{q). 



(5) 



Here L^Ly is the area of the 2DEG, = l/27rZ|j the 
number of states on one Landau level, and the form- 
factor FnIq) is 



(6) 



where Lm{x) is the Laguerre polynomial. For the case 
TV ^ 1, one can use the following asymptotic expression 
for the form-factor © 



FN{q) - MqRc) 



qRc «^ = '^ (7) 



H 



After the Hartree-Fock decouplingStt of the interaction 
term in the action |(2Jl we obtain 



S 



(8) 



T 

a 



iLUn+p-Ho- Vd^s (r) + A(r) V^" (r) 



nhiLxLy 



■ / [/(g)A(q)A(-q), (9) 
J a 



where the potential A(r) results from the perturbation of 
the uniform electron density by the charge density wave, 
and is connected with the CDW order parameter as fol- 
lows 



m = LxLyU{q)F^\q)A{(\) 



and U{q) = —n]JJHF{q) with the Hartree-Fock potential 
Unpiq) given by 



UHFiq) = Uo{q)F^{q) - 



nL 



-C/o(P)F2(p). 



(11) 



The averaging over the random potential Vdis(r) in 
Eq.(^) is straightforward and results in the following 
quartic term 



f / E ^S„(r)C„(r)V'a„(r)^f^(r) (12) 



in the action. This term can be decoupled by means 
of the Hubbard-Stratonovich transformation— with the 
Hermitian matrix field variables^SiS Q"^(7^ 



V[Q] exp 



1 

'2^ 



(13) 



where the symbol tr denotes the matrix trace over the 
Matsubara and replica indices. The measure for the func- 
tional integral over the matrix field Q is defined as: the 
integral H13() equals unity when the fermionic fields "0^ 
and "0 vanish. Note also that in Ea. (|13|l we introduce the 
matrix notations according to 



V^t(...)^ 



a,0 



/nra tuj 



(14) 



After making all these steps, the action becomes 

NrnA 1 



s 



T 2g 

+ I ip\r){iuo + p-no + \ + iQ)i}j[f), (15) 

J r 

where id is the frequency matrix ('^j^^ — ^n^nm^^ • 



C. The saddle-point in the Q field 

The Q matrix field can be naturally splitted into the 
transverse V and the longitudinal P components as fol- 
lows Q — V^^PV. The longitudinal component P has 
the block-diagonal structure in the Matsubara space, 
P"^ cx Q{nm), where Q{x) is the Heaviside step func- 
tion, and corresponds to massive modes. The transverse 
component y is a unitary rotation and describes massless 
(diffusive) modes (see Refs— for details). 

This decomposition of the variable Q into P and V is 
motivated by the saddle-point structure of the action H15|l 
at zero temperature (w„ — > 0) and in the absence of the 
potential A(r). The corresponding saddle-point solution 
has the form Qsp = V~^PspV, where the matrix Pgp is 



(p ]aP — pn Sj 



gap 



(10) with P^p obeying the equation 



^proP" -*GJf(r,r). 



(16) 



(17) 



This equation is equivalent to the self-consistent Born 
approximation equation^. The Green function G'o(r,r ) 
is determined as 



Go"(r,r 



(18) 
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Go{uJn) = [i L^n +^J■ - f^N + iP", 



(19) 



where = ujh{N + 1/2) and (j>Nk{'<^) are the eigenvalues 
and eigenfunctions of the hamihonian Tip, and k denotes 
pseudomomentum. 

In the case of a small disorder lohtq S> 1 the solution 
of equation (|17|l has the forr^ 



p:: 



2t 



P To 



m y/uJHTo 



(20) 



The fluctuations of the V field are responsible for the 
localization corrections to the conductivity (in the weak 
localization regime they correspond to the maximally 
crossed diagrams). However, in the considered case, these 
corrections are of the order of 1/A^ <C 1 and, therefore, 
can be neglected. For this reason we simply put V — 1. 

The presence of the potential A results in a shift of the 
saddle-point value (|20() due to the coupling to the fluctu- 
ations SP = P — Psp of the P field. The corresponding 
effective action for the SP field follows from Eq. after 
integrating out fermions: 



S[SP,X] = / trlnGo^ - 



^-±1 tr(P, 



5Pf 



tr In 



l + {i5P + X)G^^ 



(21) 



As a result, the thermodynamical potential can be writ- 
ten as 



T 



In j V[6P]I[6P]eyi^S[6P,\], (22) 



where, following Ref^-*, the integration measure I[5P] is 



\nI[6P] = - 



a/3 

E 



[1 - Q{nm)] SP^^SPZ.- (23) 



The quadratic in SP part of the action (|21|l together 
with the contribution ^2'^\i from the integration measure 
determine the propagator of the 6P fields (see Reii^ for 
details) 



2[l-e(mim3)] gSmirn2^'''^ gSmsmJ^'^ 



{TTpy 



l + 5^r(0;g)l-H5^o"M0;9)^ 



(24) 



r 



where the bare polarization operator 7r™(n; q) is 

TT'o{n;q) = ~nLGa{io„, + i^„)Ga{i^„^)F^{q). (25) 

D. Thermodynamical potential 

To find the expansion of the thermodynamical poten- 
tial in powers of the CDW order parameter A(q), it is 

convenient to introduce a new variable 6P = SP + iX and 
expand tr In in the action (|21|l in powers of this new field 
6P. Then the thermodynamical potential can be written 
in the form 



where 



Oo(M)=_^trlnGoi-^^trp2 



(26) 



(27) 



is the mean-field thermodynamical potential of the ho- 
mogeneous state, and 

Sn = -^\n J V[SP] exp5[^. A] (28) 

takes into account the fluctuations of the massive lon- 
gitudinal field SP and their interaction with the CDW 
order parameter (potential A). The action S[SP,X] has 
the form 



S - 5(2) [A] + S„u [SP, X] + 5(2) [SP] + J2 5^") [SP] , (29) 

n— 3 

with 



S,nt[SP,X] = -- I A(r)tr^P(r), (31) 
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and 



5^") [6P] = 



tr Yl 



<5P(r,)Go(r,r,+i), (32) 



where r^+i = Vi. Note that the terms in the action (|29|l . 
that are proportional to N^, are omitted because they 
do not contribute to in the rephca hmit Nr 0. 
Another important observation is that the propagator of 
the SP fields is the same as for the SP fields 
By using Eas. (|28|l - (|32|l we can write 



sn 



T 

A(r)A(r) - — ln(exp5i„t 



(33) 



where (• • • ) denotes the average over SP with respect 
to the action 5[(5P, 0]. This equation allows us to find 
the contributions to the thermodynamical potential 
up to any order of the CDW order parameter A(q) = 
FN{q)U{q)-'Xi(l)/L,Ly. 

In this paper we will work only with the expansion 
upto the fourth order term (the Landau expansion) . This 
implies that our consideration is valid only close to the 
transition point where the value of the order parameter 
is small and one can truncate the series (|28|l after several 
first terms. It should be mentioned, however, that we 
should avoid a direct vicinity of the phase transition (the 
critical region, for more details see Sec. llVjl where the 
fluctuations of the order parameter break the mean-field 
approach. 



1. Second-order contribution 

The second order contribution to the thermodynamical 
potential 6^1 is 

*"'^'-^i:iWAW-^(SL>„, (34) 



where (• • • )o stands for the average over SP with respect 
to the action S^^^'ISP]. We replace the average over the 
full action S[SP, 0] by the average over the quadratic part 
5^^-' [SP] only because the higher order in SP terms lead 
to the contributions that are proportional to N^, and 
therefore vanish in the replica limit iV^ 0. 

With the help of Eqs.((Tni),I21l, and we obtain 



U^{q)Gl{iO„) 
l + 5^o""(0,g) 



A(q)A(-q). (35) 



The corresponding diagram in the usual "cross tech- 
nique" is shown in Fig. ^ 




FIG. 1: Second-order contribution to the thermodynamic 
potential. Solid line denotes electron Green function, dashes 
are impurity lines and vertexes are A(r) 




FIG. 2: The third-order contribution to the thermodynamical 
potential. 



2. Third-order contribution 

The contribution of the third power of the CDW order 
parameter to the thermodynamical potential Jfi'^' can 
be written as 



T v(c) 



3!iV, 



/q3 y^) ^ 

\'-'mt/ xp 



T 



3,^ \SLS^^'^)^ , (36) 



where the superscript (c) indicates that only connected 
diagrams are taken into account. Here we omit again 
the terms that vanish in the replica limit Nr ^ 0. Af- 
ter performing the averaging over SP with the help of 
Eqs.|(in|l,(|2iI),lElll, and we obtain 



UJ„ j = l 



U{q,)A{q,)Go{uJn) 
q, l+57ro"(0,gj) 



X (5(qi+q2 + q3)exp^((7j'g|-gfgf). (37) 

The contribution Sfl'^^^ corresponds to the diagram in 
Fig. 12 



3. Forth-order contribution 



The forth order contribution (517 is 



T ,(c) 



T 



c4 



1 



5(4) + ±(5(3))2 \ ^ (38^ 

2 J / Q 



6 




FIG. 3: The forth-order contribution to the thermodynamical 
potential. 



where is the free energy of the normal (homogeneous) 
state, Nf. the total number of electrons, and //q the 
chemical potentials of the CDW state and the normal 
state respectively. 

In order to find the free energy of the CDW state to 
the forth order in the CDW order parameter we expand 
^o(mo) around the point to the second order in /i — /iq. 
This results in 



where again only terms which is proportional to iV^ are 
kept. By using Eqs.|(ini),lE3,lElll, and we find 



T = T, + n{|,)-n,{^l)^\{^,-^,of^. (4i) 

2 o^/xo 



- (2-)v-En 



j7(g,-)A(q,)GoK) 



X (5(qi + q2 + q3 + q4) 



i-<?<"(o,|qi + q2 



i + 5<"(o,|qi + q2|) 

X exp ^ (qfql^ - qfq^'j exp ^ (q^qf - q^qfj (39) 

In the usual "cross technique" the contribution Sfl^^'' cor- 
responds to the diagram shown in Fig. |2| 



The difference /i — /ig of the chemical potentials in the 
CDW and the normal states is 



M - Mo 



dsn fdNe\~'^ 



and from Ea. H26(l we obtain 



(42) 



E. Free energy 

The free energy of the CDW state can be written in 
the form 



^^^, + n. + in + i(^l7^V'. (43) 



2 \ J \ dfi 



J- = J-Q + fi(/^) — rio(Mo) + (p ~ Mo)-^e, (40) With the expression ^ for this gives 



A ■ 



(50 



nL{LxLy) 



T 



UHq)GlM 
, [l+5^S'"(0,g)]^ 



■A(q)A(-q) 



(44) 



r 



F. Free energy of the triangular CDW state 

The CDW order parameter for the triangular lattice 
symmetry (bubble phase) can be written in the formi 



expansion are as follows 

ToiQ) 



a2 — 3 



1 - 



E 



(47) 



A(q) 



(2^ 

LxLy 



A(g)^[<5( 



q-Q,) + <5(q+Q,) , (45) 



where the vectors have the angle 27r/3 between each 
other and obey the condition Qj^ + Q2 + Q3 = 0. 

By using Eqs. 11231,(123,10 and we obtain the 

following expression for the free energy of the triangular 
CDW state 



. LxLy 

2^ 



a2A^ + asA^ + a4A^ 



(46) 



where 



1 



1 



2 4ttTt 2ttT 
andTo(Q) = C/(g)/4, 



as = i( 



7r3T2 



Here the three coefficients ai , 02 , and 03 of the Landau and 



(48) 



(49) 



7 



247f(g)/l>^ 



+ 3 



3A.(0) 



1 + COS ) f DniQ) + DniVSQ) 



-M2Q) 



(50) 



with 



(51) 



G. Free energy of the unidirectional CDW state 

The CDW order parameter of the unidirectional state 
(stripe phase) isSiS 



A(q) 



(2^)^ 
LxLy 



A(Q) (5(q-Q)+(5(q-Q) 



(52) 



where the vector Q is oriented along the spontaneously 
chosen direction, and from Eqs. (123,(1237 (123) and (g3, 



the free energy of the unidirectional CDW state reads 



(53) 



Here the coefficients 62 and 64 of the Landau expansion 



h °2 
02 = -7- 



and 



E 



1 



ToiQ) 



(54) 



64 = 



D^{Q) + -D^{2Q) 



E 



[EC 



(55) 



We note that in the limit 1/r ^ evaluation of fre- 
quency sums in expressions for free energy of the trian- 
gular and unidirectional CDW states 103- (1^3 leads to 
the results obtained in Refsi^ in the clean case. 



tion 



TT 



(56) 



III. MEAN-FIELD PHASE DIAGRAM 

A. Instability line 

The vanishing of the coefficient in front of the 
quadratic term in the Landau expansion of the free en- 
ergy signals about the instability of the normal state to- 
wards the formation of the CDW. This instability corre- 
sponds to the second order phase transition from the ho- 
mogeneous state to the CDW state. As usual, the specific 
parameters of the forming CDW state are determined by 
the high order terms in the Landau expansion. 

From Eqs. H47|) and H54|l we obtain the following equa- 



for the instability line. The solution T{Q) of this equa- 
tion depends on the modulus Q of the vector that charac- 
terizes the CDW state. The temperature T2 of the second 
order phase transition corresponds to the maximal value 
of r(Q): 



T2 =max T(Q) 



(57) 



and the corresponding value Qo, T2 = T{Qo), determines 
the period of the CDW state. The Hartree-Fock poten- 
tial 1)11(1 has minima at those vectors Qk for which the 
form-factor F]s[{Qk) vanishes. In the clean case this cor- 
responds to Qo = minQfc = tq/Rc, where tq « 2.4 is the 
first zero of the Bessel function of the first kind—. It can 
be seen from Eq. H56|l that a weak disorder does not shift 
the vector Qo (see Appendix) . Thus the equation for the 
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temperature of the second order phase transition into the 
CDW state reads 



T 2 „ , /I 1 
To TT^ ^ \2 AttTt 



+ i 



2t:T 



(58) 



where ip (z) is the derivative of digamma function, 3? 
the real part, and Tq = Tq{Qo) the temperature of the 
transition in the clean case. 

Ea. H58|) contains the chemical potential /i that, to- 
gether with the temperature T and the broadening of 
Landau levels 1/2t, determines the filling factor i/^ = 
h' — 2N of the partially filled highest Landau level. How- 
ever, in order to find this relation one needs to know 
the density of states in the system. This question about 
the density of states is a very subtleSi and beyond the 
scope of the present paper. For this we use the chemical 
potential /i rather than filling factor v^. 

Ea. H58|) can be solved analytically in the two extreme 
cases; when temperature T is closed to the temperature 
To of instability in the absence of disorder, and when the 
temperature T is close to zero. 

In the first case, the broadening of the Landau level 
1/2t and the chemical potential /x are small compared to 
the temperature Tq of the instability in the clean case, 
and, therefore, the leading order expansion in powers of 
1/Tqt and /i/To is legitimate. It appears that the pres- 
ence of disorder decreases the temperature of instability 
linearly: 



T 



7<(3) 



4T2 



In the opposite case T 
and Ea. l(5^ reduces to 



— ,M«2^r. (59) 
0, one has 1/2t, ^ > 27rT, 



1 



8Tqt 1 + 4^^^ 



(60) 



We see from Ea. (|60f) that at zero temperature the second 
order phase transition can occur only when the broaden- 
ing of the Landau level is smaller than some critical value, 
1/2t < l/2rc = iTo/iT. 

For other cases Eq. H58|) can be solved numerically, and 
the corresponding instability (spinodal) line is shown in 
Fig.ll 



B. Half-filled Landau level (un = 1/2) 



I- 

*^ 

m 

L. 

o 
a 



J I I I I I L. 





= 0.0 


jt/8T^T 


= 0.1 


mi^x 


= 0.5 


mi^x 


= 0.9 




— I — •- 

0.6 0.8 1.0 1.2 

chemical potential, \il 2iil 



1.4 



1.6 



FIG. 4: The spinodal lines obtained from Eg. (1581 are shown 
for different values of dimensionless parameter tt/STot. 
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Landau level index = 2, 3, 4, 5 increases from left to right. 



We now consider the case of the half-filled A''th Lan- 
dau level {vn = 1/2), that is related to the recent experi- 
ments^. In this case the chemical potential is zero, /i = 0, 
provided the density of states is symmetric around the 
center of the TVth Landau level. As follows from Ea. (|58|l . 
the temperature of the second order phase transition for 
this case can be found from the equation 

T 2/1 1 \ 



where ^(2, z) = J2m=o(^ ~^ ^)^^ is the generalized Rie- 
mann zeta function. The analytical solutions of this 
equation in the cases of high and low temperature can 
be obtained from Eas. H59|l and 1)60(1 by putting n to zero. 
The entire behavior of the spinodal line, obtained numer- 
ically from Ea. l|61|) . is shown in Fig.|Sl 

We mention that at = 1/2 the coefficient van- 
ishes due to the particle-hole symmetry. It means that 
the transition from the normal state into the CDW state 



9 



1.00 



0.95 



I- 

co 



o 
■o 



S 0.90 



0.85 



J I I I I I I 1 I I L 





Uniform 
^ phase 

s N. 


Unidirectional 
CDW phase 

— 1 — 1 — 1 — 1 — 1 — 1 — 1 


\ N. 
S 

\ \ 
\ \ 
\ \ 
\ \ 
\ \ 
\ \ 
\ \ 

Triangular \ \ 
CDW phase \ 

\ 

\ 

1 1 1 1 1 1 1 1 1 1 |V 



0.00 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 

Chemical potential, |a x 

FIG. 6: Phase diagram at zero temperature near vn ~ 1/2. 
The sohd line is obtained from Eg. 1)72^ . the dashes are the 
spinodal line and the dots are obtained from Ea. l74t . 



is of the second order for both cases of unidirectional and 
triangular lattice symmetry. Therefore, to find the struc- 
ture of the CDW state, one has to take into account the 
fourth order terms in the the Landau expansion. In the 
vicinity of the spinodal line, it follows from Eqs.ljSOJ and 
with /i = that 



04 : 

and 



12T§ 



-7C(4, u) + 12$o(w) + 6$2(u) + 8$^(u) 



2T3 



4y3 



-3C(4,m) + 4$o(") + 2$2(m) 



(62) 
(63) 



where we introduce the new variable u = 1/2 + I/AttTt 
and the new function 



1 



1 



C(2,i + z 



(64) 



^■ip - + z + i zjQ{arQ) 



zJoiaro) 



with 5 being the imaginary part. With these expressions 
we minimized J^*^" with respect to the order parameter A 
and found that the unidirectional CDW state has lower 
free energy. 



C. Phase diagram at zero temperature 

In this section we analyze the zero temperature phase 
diagram in the case where the Landau level broadeding 



is close to its critical value 1/2tc = ATq/tt. Under these 
conditions, the CDW order parameter A is small, and 
one can use the Landau expansions (|46l) and l|53|l at zero 
temperature. The coefficients of these expansions are 



aa = 3 1 



and 



(65) 



«4 = 3^'(^^j H3ifST,QR,) 
for the triangular CDW state, and 



(66) 



62 = 



a2 



64 = 



Hi{nT,QRa) (67) 



3 ' 2 V TT 

for the unidirectional CDW state. Here wc introduce four 
fimctions Hi{z) as 

1 



iJi(z) = 



1+4^2 



H2{z) = izHiiz), 



(68) 



H3iz,r) 

and 

Hi{z,r) 
where 

Ra{z,r) = 



1 108z^ - 5 , , , , 

2 (l + 4z2)3 +3-«o(^,r) + 2i?i(z,r) 

lR2{z,r) + 2R^iz,r), (69) 



28z" 



1 



(l + 4z2)3 

2iJi(z) 



+ 2Ra{z,r)+R2iz,r), (70) 

2Jo(ar) 
l + Az^-j;^{ar)' 



arctan 



These expressions result in the following equation on 
the line of the first order transition from the uniform to 
the triangular CDW state 



8Tnr 



(72) 



9i/3(A*T, QRc)' 

As before, the maximum of the solutions l/r(Q) of 
Eq.lO with respect to Q should be found. It appears 
that the maximum is not exactly at Q = Qq as in the 
clean case^, but at some shifted value Qo + SQ with the 
shift SQ = -0.02(^r)2i?ji for small fir < 1. The ex- 
istence of the shift is a feature of the disordered case. 
Below in the limit //t <C 1, we will neglect this shift. In 
this case Eg . (1721) can be written as 



= 1 - 2.94(/ir)2 



(73) 



STor 

By comparing the free energies of the triangular and 
the unidirectional CDW states, we can find the line of 
the first order transition between them 



8T0T 



Hi 



Hi[2H3 + Hi][dHi + ^Hl + 2i/3i?4] 



2^3(2^3 -3i/4)2 



(74) 
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For the case /ir ^ 1 Eg . (|74() can be simplified as 
= 1- 18.44(^r)2. 



8ToT 



(75) 



For other values of /ir Eas. l|72|l and Ea. l74|) were solved 
numerically, and the results are shown in Fig. El When 
the parameter S/ttTqt decreases at a fixed value of the 
chemical potential, the CDW order parameter grows, and 
hence, we go beyond the applicability of to the Landau 
expansion. 



IV. WEAK CRYSTALLIZATION 
CORRECTIONS 

The CDW order parameter A(r) introduced in Eq.© 
can be thought of as a saddle-point solution for the 
plasmon field that appears in the Hubbard-Stratonovich 
transformation of the electron-electron interaction in the 
action (j^J. The Landau expansions (|46|l and H53(l for the 
free energy of the CDW states were derived under the 
assumption that one can neglect the fluctuations of the 
CDW order parameter. This is legitimate for iV ^ 1 
and not very close to the transition (outside the critical 
region). However, when one approaches the instability 
line, the fluctuations of the CDW order parameter in- 
crease. To analyse the effects of the order parameter 
fluctuation, we introduce, following the original ideas of 
Brazovski^^, the fluctuations of the CDW order parame- 
ter A(r) A(r) + S{r) in the Landau expansion of the 
free energy and average over the fluctuations 6{r). We 
present below the results of the corresponding analysis 
only for the most interesting case of the half-filled Lan- 
dau level. 

We find that the transition from the uniform to the 
unidirectional CDW state becomes of the first order, and 
takes place at the lower temperature that can be found 
from the following equation (see Eq. (|61l for comparison ) 



T 2/1 
^ = — C 2,- 



1 



( 1 



2 ' 47rrT ) " \ AtiTt 
Here function g{z) is defined as 

where we introduce the following three functions 



iV"^/^ (76) 



"37rro' 


2/3 




2/3 


■2Ao(z) + 


A2(^)l 


1/3 


16 




[ /(^) \ 






A2(^)J 





/(-) = Z2 



<(4,i + .)+2a>. (4,1 



(77) 
(78) 

(79) 



The constants [3i are given by 



ToiQo) 



2.58 , p2 = (Jo'(Qo)) « 0.27, (80) 



and function <I>a is defined by Ea. l|^ . 

We mention that the function 7r^(7(z)/2C(2, 1/2 + z) 
decreases monotonically from the value 0.35 at z = 
to zero at z — > cxD. Therefore, we obtain the following 
inequality for the shift 5T of the mean-field transition 
temperature T 



ST , 
— < 3 , 
T - Vl6 



2/3 



iV>l (81) 



(the equality corresponds to the clean case. 

The appearance of a noninteger powers in Ea. H76|l re- 
sults from the fact that the momentum dependence of the 
correlation function for the order parameter fluctuations 
contains {Q — Qo)^ rather than (see Refi^). 

Ea. l|76|l was derived under the assumption that the 
main contribution in the momentum space comes from 
the region Q « Qq. This assumption is justified under 
the following condition^* 



inTr 



1 



(82) 



The combination of functions in the left hand side of 
inequality H82|) decreases monotonically from 0.023 to 
while z increases from zero to infinity and, hence, the 
condition (|82|l is hold. 

According to Ea. H81|l . the fiuctuations reduce the 
transition temperature by the amount of the order of 
^-2/3 ^ Y slTlA, therefore, in the considered case of the 
weak magnetic field {N ^ 1) their effects can be ne- 
glected. These results indicate that the critical region for 
the considered transition is indeed small, and the mean- 
field approach gives a good approximation for N ^ \. 



V. DISCUSSIONS 

A. Comparison with experimental results 

Now we discuss the possible applications of our theory 
to the recent experiments. Although our mean-field the- 
ory was derived for the case of a large number of the occu- 
pied Landau level N 1, and neglects corrections of the 
order of 1/A^, while experimentally one has N = 2,3,4, 
we however expect that Ea. (|58|l gives a good estimation 
for the temperature of the transition from the uniform 
to the CDW state, even for TV = 2,3,4. We have com- 
plementary assurance that it can really be the case be- 
cause Ea. (|58|l can be obtained without introducing the 
CDW order parameter and considering the mean-field 
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theory but as the equation that determines the temper- 
ature T(Q) at which the two-particle vertex function at 
wave vector Q diverges^S. 

We restrict ourselves by discussion of the experiments 
without an in-plane magnetic fieldSiiS. The theory for 
the half-filled highest Landau level contains two physical 
parameters: the temperature Tq and the broadening of 
the Landau level 1/2t. The estimate of the temperature 
of instability Tq in the absence of disorder is more sub- 
tle. The theory can provide an estimate for U{Qq) and, 
correspondingly, for the value of Tq only for a weak mag- 
netic field. In this case the value of U{Qq) can be found 
analytically^ with the following result for Tq 

To = ^LUH , NTs » 1 , < I, (83) 

where the numerical parameter a only slightly depends 
on Vs in the range 0.1 < < 1, and approximately 
a ~ 0.03. The broadening of the Landau level 1/2t can 
be estimated from the mobility at zero magnetic field. 
However, the results obtained in Sec. IIII Bl impose the 
restriction on the value of the sample mobility at zero 
magnetic field. In order to observe the CDW states at 
the partially filled Landau level with index N the mo- 
bility at zero magnetic field should satisfy the condition 
/i > eN I2c?ne where rie denotes the electron density of 
2DEG. For typical values of the electron density of 2DEG 
we can obtain the following estimate n > N ■ lO^cm^ /Vs. 

In the experiments of Lilly et aliS, the samples were rel- 
atively clean (mobilities exceed 9 • 10^ cm^/Vs), and the 
anisotropy was observed for N — 2 below T — 150 niK. 
For the temperature T = 25 mK, the anisotropy in the 
resistance was found to disappear with the decreasing the 
magnetic field for iV = 4. We plot experimental points in 
Fig. [31 under the assumption that Eq. H83|l remains valid 
even for N = 2,3,4, and 5. As it can be seen from this 
Figure, for a constant magnetic field the highest tempera- 
ture at which anisotropy appears, corresponds to TV = 2. 
As the magnetic field decreases at some constant temper- 
ature, the anisotropy disappears at some value of TV due 
to the disorder induced transition from the CDW state 
to the uniform. It should be mentioned that without the 
effects of disorder the anisotropy in the resistance should 
remains up to iV = 12 at T = 2bmK. Therefore, the dis- 
order plays an important role even in the high mobility 
samples. 

B. Comparison with numerical results 

The problem of the formation of the CDW state on 
the second Landau level with i^at = 1/2 at zero temper- 
ature in the presence of a quenched disorder was studied 
numerically in Refi^S. The system of 12 electrons inter- 
acting via the Coulomb interaction U{q) = 2Tre'^ /q in the 
presence of the quenched disorder was projected on the 
second Landau level {N = 2). The effects of interactions 
with electrons on the other Landau level was not taken 



into account. The system was diagonalized numerically. 
It was found that the CDW state transforms into a uni- 
form liquid state as the dimensionless disorder strength 
i-^H\/nL /27rpTo exceeds 0.12. 

In order to be able to compare the results of the pre- 
sented above mean-field theory with the numerical re- 
sults, we perform the evaluation of the temperature Tq 
in the case for which the numerical results were obtained 
(instead screened interaction (|3J) we use U{q) = 27re^/g). 
Under this curcumstances our theory gives the value 0.14. 

The small discrepancy may be attributed to two fac- 
tors: on the one hand, the finite number of electrons 
in numerical calculations and, on the one hand, unsufR- 
cience of the Hartree-Fock approximation for the prob- 
lem with the Coulomb interaction U{q) = 2'Ke^ jq. In 
the later case one should take into account the diagrams 
beyond the Hartree-Fock theory. Nevertheless, the com- 
parison demonstrates that such corrections are small. 

We emphasize that our theory which takes into account 
the screening of electron-electron interaction by electrons 
on the other Landau levels gives much smaller value 0.038 
of the dimensionless disorder strength for the transition 
from the uniform liquid state to the CDW state. 

VI. CONCLUSIONS 

For the system of a two-dimensional interacting elec- 
trons in the presence of a weak disorder and a weak mag- 
netic field, we investigated the effect of disorder on the 
existence of the mean-field CDW states in the frame- 
work of the Hartree-Fock approximation. In the consid- 
ered case of large filling factors ^ 1, we obtained that 
the mean-field CDW instability exists if the disorder is 
rather weak, 1/r < STq/tt. We found that at half- filling 
the unidirectional CDW state appears, and the presence 
of disorder does not change the vector of the CDW. Near 
half-filling, the unidirectional CDW state is energetically 
more favorable than the triangular one. We obtained that 
the weak crystallization corrections to the mean-field re- 
sult are of the order of (1/z/)^/'^ ^ 1 and thus can be 
neglected. We discussed the applications of our theory 
to the recent experimental and numerical results. 

Acknowledgments 

LB. is grateful to M.V. Feigel'man and M.A. Skvortsov 
for stimulating discussions. This research was sup- 
ported by Forschungszcntrum Jiilich (Landau Scholar- 
ship), Russian Foundation for Basic Research (RFBR), 
and Deutsche Forschungsgemeinschaft (DFG). 

APPENDIX A: INSTABILITY VECTOR 

In this appendix we prove that the weak disorder does 
not change the vector at which the instability towards the 
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formation of the CDW state grows. Let us consider the 
solution T + ST of Eq.lO for the vector Q ^ Qo + SQ, 
where SQ <C Qo- We wiU now show that the shift ST 
is always negative, and hence, the maximal instability 
temperature corresponds to the vector Q — Qo, as it is 
in the clean case. 

For a small deviation SQ we can write 

ToiQ) = ro(l - Pii5QRcf) , J^iQRc) = p2{5QR,f. 

(Al) 

The shift 5T results in the substitution 



5T 



.5TjJ_ 



in Ea. H56() . and we obtain 



— = -l^dQHc) ) — ■ =r, 

T 52 - 2{zg3 + yga) 



(A2) 



(A3) 



where z = I/AttTt and y = Here we introduce 

the four functions ga{z,y) and ga{z,y) 



n=0 



n=0 



It can be easily seen that the rhs of Eq. (|A3p is negative 
for all possible values of z and y. 



* Also at Kurchatov Institute, Kurchatov Square 1, 123182 
Moscow, Russia 

^ For a review, see The Quantum Hall effect, ed. by R.E. 

Prange and S.M. Girvin (Springer- Verlag, Berlin, 1987) 
^ A.H. MacDonald and S.M. Girvin, Phys. Rev. B 33, 4009 

(1986) 

^ L. Belkhir and J.K. Jain, Solid State Commun. 94, 107 
(1995); R. Morf and N. d'Ambrumenil, Phys. Rev. Lett. 
74, 5116 (1995) 

* A.P. Smith, A.H. MacDonald and G. Gumbs, Phys. Rev. 
B 45, 8829 (1992) 

^ I.L. Aleiner and L.I. Glazman, Phys. Rev. B 52, 11296 
(1995) 

® A.A. Koulakov, M.M. Fogler, and B.I. Shklovskii, Phys. 

Rev. Lett. 76, 499 (1996), Phys. Rev. B. 54, 1853 (1996) 

H. Fukuyama, P.M. Platzmann, and P.W. Anderson, Phys. 

Rev. B 19, 5211 (1979) 
^ R. Moessner and J.T. Chalker, Phys. Rev. B 54, 5006 

(1996) 

^ M.P. Lilly, K.B. Cooper, J. P. Eisenstein, L.N. Pfeiffer, and 

K.W. West, Phys. Rev. Lett. 82, 394 (1999) 
1° R.R. Du, D.C. Tsui, H.L. Stormer, L.N. Pfeiffer, and K.W. 
West, Solid State Commun. 109, 389 (1999) 
For a review, see M.M. Fogler, arXiv: cond-mat/ 0111001 
and references therein. 

A. Lopatnikova, S.H. Simon, B.I. Halperin, and X.-G. Wen, 
Phys. Rev. B 64, 155301 (2002) 

D.G. Barci, E. Fradkin, S.A. Kivelson, and V. Oganesyan, 
Phys. Rev. B 65, 245319 (2002); D.G. Barci and E. Frad- 
kin, Phys. Rev. B 65, 245320 (2002) 
^'^ S. Scheidl and F. von Oppen, Europhys. Lett. 55, 260 
(2001) 

W. Pan, R.R. Du, H.L. Stormer, D.C. Tsui, L.N. Pfeif- 
fer, K.W. Baldwin, and K.W. West, Phys. Rev. Lett. 83, 



820, (1999); M.P. Lilly, K.B. Cooper, J.P. Eisenstein, L.N. 
Pfeiffer, and K.W. West, Phys. Rev. Lett. 83, 824, (1999) 
A.A. Abrikosov and L.P. Gor'kov, Zh. Eksp. Teor. Fiz. 39, 
1781 (1960) 

A.I. Larkin, Pis'ma Zh. Eksp. Teor. Fiz. 2, 205 (1965) 
I.S. Burmistrov, Zh. Eksp. Teor. Fiz. 122, 150 (2002) [ 
JETP 95, 132 (2002) ], arXiv: cond-mat/0203022 
X.-G. Wu and S.L. Sonhdi, Phys. Rev. B 51, 14725 (1995) 
For details of the HF decoupling see e.g. Reffi and Refj^. 
R.L. Stratonovich, Dokl. Akad. Nauk SSSR 2, 1097 (1957) 
[Sov. Phys. Doklady 2, 416 (1957)]; J. Hubbard, Phys. Rev. 
Lett. 3, 77 (1959) 

K.B. Efetov, A.I. Larkin, D.E. Khemel'nitzkii, Zh. Eksp. 
Teor. Fiz. 79, 1120 (1980)[Sov. Phys. JETP , (1980) 
A.M. Finkel'stein, Pis'ma Zh. Eksp. Teor. Fiz. 37, 436 
(1983) [ JETP Lett. 37, 517 (1983) ]; Zh. Eksp. Teor. Fiz. 
84, 168 (1983) [ Sov. Phys. JETP 57, 97 (1983) ]; Zh. 
Eksp. Teor. Fiz. 86, 367 (1984) [ Sov. Phys. JETP 459, 
212 (1984) ] 

A.M.M. Pruisken, Nucl. Phys. B 235, 277 (1984) 
A.M.M. Pruisken, M.A. Baranov, and B. Skoric, Phys. 
Rev. B 60, 16807 (1999) 

T. Ando and Y. Uemura, J. Phys. Soc. Japan 36, 959 
(1974); tbid 36, 1521 (1974); T. Ando, ibid 37, 1233 (1974) 
For a review, see I.V. Kukushkin, S.V. Meshkov, and V.B. 
Timofeev, Usp. Fiz. Nauk 155, 219 (1988) [Sov. Phys. Usp. 
31, 511, (1988)] 

S.A. Brazovski, Zh. Eksp. Teor. Fiz. 68, 175 (1975) [Sov. 
Phys. JETP, 41, 85 (1975)] 

for the clean case see Ref.** and for the weakly disordered 
case I.S. Burmistrov and M.A. Baranov unpublished 
D.N. Sheng, Z. Wang, and B. Friedman, Phys. Rev. B 64, 
195316 (2002) 



